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Abstract. The distribution of gravitational torques and bar strengths in the local Universe is derived from a 
detailed study of 163 galaxies observed in the near-infrared. The results are compared with numerical models 
for spiral galaxy evolution. It is found that the observed distribution of torques can be accounted for only with 
external accretion of gas onto spiral disks. Accretion is responsible for bar renewal - after the dissolution of 
primordial bars - as well as the maintenance of spiral structures. Models of isolated, non-accreting galaxies are 
ruled out. Moderate accretion rates do not explain the observational results: it is shown that galactic disks should 
double their mass in less than the Hubble time. The best fit is obtained if spiral galaxies are open systems, still 
forming today by continuous gas accretion, doubling their mass every 10 billion years. 
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^ ' 1. Introduction 



Bars are a major perturbation of the gravitational po- 
tential and a highly efficient engine for the evolution of 
morphological and chemical properties. In purely stellar 
disks, bars are robust, long-lived structures (Combes & 
Sanders 1981). In reality, however, spiral galaxies contain 
gas which provokes the dissolution of bars (Bournaud & 
Combes 2002, hereafter BC). The bar itself initiates an 
important radial gas inflow, which destroys barred struc- 
ture (Pfenniger & Norman 1990). 

Since bars in gaseous disks are dissolved in a few Gyrs, 
the ubiquity of bars in the local Universe suggests that 
bars are also renewed. BC painted a scenario wherein con- 
tinual gas accretion might be a crucial evolutionary factor. 
This gas settles in the disk, enhancing its self gravity, re- 
ducing the influence of the stabilizing central mass concen- 
tration, and a second bar develops. In this Letter, we test 
this hypothesis, by comparing the observed distribution 
of gravitational torques (in a well-defined galaxy sample) 
with the distribution predicted by numerical simulations 
incorporating gas accretion. 



2. Observational and numerical techniques 

2.1. Gravitational torques 

The maximal torque over the entire disk is defined by 



77' max 



1 f d<S>(RM) ^ 

R\ ae I 



d<!>o(R) 
dR 



(1) 



where F™'''' 



Send offprint requests 
blockOcam. wits . ac . za 



to: 



D. 



Block e-mail: 



(i?) represents the maximum amplitude of the 
tangential force and Fq(R) is the mean axisymmetric ra- 
dial force, inferred from the axisymmetric component of 
the gravitational potential, $o- The potential is inferred 
from the visible mass only, in both simulations and ob- 
servations. We then measure Qh — ina.xji{QT{R)). Qh is 
simply the bar strength in barred galaxies, or the maximal 
arms torque in unbarred or nearly unbarred galaxies. 

2.2. Observations 

To derive the distribution of Qb, we use near- infrared im- 
ages from the Ohio State University Bright Galaxy Survey 
(OSUBGS, Eskridge et al. 2002). The full sample consists 
of 198 spiral galaxies having total apparent blue magni- 
tude Bt < 12 and blue isophotal diameter D25 < 6'. The 
technique for deriving Qb from near-IR images is described 
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Type 


bulge-to-disk 


halo-to-disk 


abundance 




mass ratio 


mass ratio 


(OSU sample) 


Sa 


0.9 


0.3 


17% 


Sb 


0.4 


0.5 


25% 


Sc 


0.2 


0.7 


45% 


Sd 


0.1 


0.9 


13% 



Table 1. Bulge, disk, and halo parameters used in nu- 
merical simulations. The halo mass that is accounted for 
is the mass of dark matter inside the disk radius. 
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by Buta & Block (2001; see also Block et al. 2001), who 
used the method of Quillen, Frogel & Gonzalez (1994) to 
transform deprojected, cleaned images into potentials. 

Of the 198 OSUBGS galaxies, only 163 are used here 
for studies of gravitational torques, for reasons detailed 



0.2 



0.4 



0.6 



0.8 



in Sect. 4.1, A constant mass-to-light ratio is assumed, 
i.e, the dark halo is ignored. The disk scale-height and 
scale-length are assumed to be related by — h^/12. 
Also, bulges are assumed to be as flat as the disk; no 
decompositions to account for the likely rounder shapes 
of some bulges have been carried out for this study. Such 
refinements will be considered in a separate study (Buta 
et al. 2002). 

2.3. Numerical simulations 

We employ a numerical code to study the evolution of bars 
and spiral arms in galaxies which do, or do not, accrete 
gas. The simulations include gas dynamics, star formation, 
and stellar mass-loss. They model a disk, a bulge, and a 
dark halo. For a complete description, see BC. 

We simulate the evolution of four model galaxies, as- 
sumed to represent types Sa, Sb, Sc and Sd, with param- 
eters given in Table |^. The proportion of each morpho- 
logical type in the observed sample is fully accounted for 
when we analyze the results. We compute the evolution of 
each model over 15 Gyr. We calculate the mean value of 
Qb at periods of 100 Myr. 

3. Results 

3.1. Distribution of gravitational torques 

A histogram showing the number of galaxies as a function 
of their maximal torque is presented in Fig. |^. The most 
striking feature of the histogram is the depression in the 
number of galaxies with very low torques (Qb ~ 0). That 
the drop at Qb is not due to selection effects, but is 
indeed robust, is discussed in Sect. 11. 

In Fig. |l|, one clearly sees that the shape of the ob- 
served histogram is not reproduced without gas accre- 
tion. Only large accretion rates reproduce the observed 
histogram (both the drop at Qb and the number of 
strong bars with Qb > 0.5), when the disk mass doubles 
in 13, 10 or 8 Gyrs. 

We also compute the deviation between the observed 
and simulated histograms. If Ni is the number of galaxies 



Fig. 1. Histogram of gravitational torques. Solid: observa- 
tions - Dashed: simulations with gas accretion that double 
the disk mass in 10 Gyrs - Dotted: simulations without 
gas accretion. 



Accretion rate 


Qb 


Qr 


Q2 


Q2 




no accretion 


1.00 


1.00 


1.00 


1.00 


1.00 


Mo/30 


0.90 


0.85 


0.86 


0.92 


0.89 


Mo/13 


0.37 


0.45 


0.51 


0.39 


0.32 


Mo/ 10 


0.34 


0.44 


0.48 


0.34 


0.26 


Mo/8 


0.40 


0.88 


0.52 


0.36 


0.28 



Table 2. Values of a for different accretion rates (in 
Gyr^^), and several gravitational torques and bar strength 
indicators. Mq is the initial disk mass. 



observed in each torque class i,and N[ the corresponding 
value from simulations, we define = J^ii^i ~ ■ We 
normalize cr to 1 for non accreting models. Values of a are 
listed in Table The best fit to the observed histogram is 
obtained when galaxies double their mass every 10 Gyrs, 
while low accretion rates in which the disk mass doubles 
in 30 Gyrs, are ruled out. 



3.2. Interpretation: temporal bar evolution 

The presence of gas in galactic disks is responsible for the 
destruction of bars in no more than 5 Gyrs. As elucidated 
by BC, gas is also responsible for bar renewal, when it 
is accreted from outside the disk. Gas accretion radically 
changes the temporal evolution of the bar strength: see 
for example Fig. 4 in BC. With accretion, once a bar is 
dissolved, the disk can become unstable again and a new 
bar may form. The disk rarely resides in axisymmetric 
states: even between bar episodes, spiral arms maintain a 
significant torque, for accretion also rejuvenates the spiral 
structure. Without accretion, the disks spend half of their 
lifetimes in nearly axisymmetric states with Qb < 0.05: 
both bars and arms disappear. The number of galaxies 
in each class of the histogram in Fig. |l| is interpreted as 
the time fraction galaxies spend in each class during their 
lives. 
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4. Discussion 

4.1. Determination of gravitational torques 

In deriving gravitational torques, we assume that each disk 
has a constant scale-height. What we measure is thus dif- 
ferent from the actual torque in nature. However, we com- 
pute Qb exactly in the same manner for the simulations. 
We adopt the same set of assumptions, which allows a fair 
comparison between observed and simulated distributions 
to be made. Even if the value of Qh is different from the 
actual maximal torque, the same error will follow through 
in our simulations. 

Obviously, Qb will differ from the actual torque also 
because of the dark matter potential. However, the max- 
imal torque is obtained well inside the disk, where the 
dark matter contribution is negligible. Moreover, Qb is 
also computed in simulations ignoring dark matter. 

A differentiation between bar and spiral arm torques is 
not needed. Spiral structure is also present in simulations, 
so knowing whether Qb is related to a bar or to a spiral 
structure is not required. 

As mentioned above, only 163 galaxies from 198 in 
the OSU sample were retained. Our selection criteria are 
delineated as follows: 

— galaxies inclined more than 70 degrees cannot be prop- 
erly deprojected (deprojection is discussed later) 

— members of obvious interacting systems are not stud- 
ied 

It is unlikely that these criteria are responsible for the 
depression in the number of galaxies with low torques: in 
particular, unbarred galaxies were not given preferential 
rejection. 

Thus far, we have itemized several approximations in 
Qb) but the same assumptions have followed through in 
the simulations, so that the comparison remains fair. The 
most questionable point is in fact deprojection, for only 
observations require such a treatment. An axisymmetric 
disk could appear non-axisymmetric when the deprojec- 
tion is not correct: here some approximations may be in- 
troduced in observations, but not in simulations. First, 
observational deprojection assumes that outer disks are 
circular, which simulations suggest is unlikely. To esti- 
mate the resulting uncertainty, we experimented by ran- 
domly orienting the simulated images, and deprojecting 
them exactly as was done for observations. If we denote 
by AQb the change introduced on each bar torque value, 
we find from simulations that \/ < AQb^ > = .022, and 
< AQb >= .007. A second problem is related to the ob- 
servational uncertainty of the disk axis ratio b/a. We use 
data from LED A to determine the uncertainty on 6/a as a 
function of b/a. We then convert the result into an uncer- 
tainty upon Qb as a function of the estimated inclination 
i of the disk, and add to this the first uncertainty on Qb, 
related to the outer disk shape; values are given in Table ||. 
Given that the bin of our histogram is 0.05, the inclina- 
tion cutoff at 70 degrees is justified, and the results of 



i 


AQb 


i 


AQb 


i 


AQb 


20 


0.022 


50 


0.025 


70 


0.039 


35 


0.023 


60 


0.028 


80 


0.052 



Table 3. Deprojection uncertainties on gravitational 
torques AQb, as a function of the disk inclination i. 



deprojection are dependable. Another deprojection issue 
is related to the bar orientation relative to the line of sight: 
this may affect strongly barred galaxies, but will not affect 
the crucial number of galaxies with low torques. 

4.2. Ottier measures of torques and bar strengths 

It is interesting to examine results obtained with physical 
indicators other than Qb. As mentioned above, Qb is re- 
lated either to a bar or to arms, for we select the maximal 
torque. It may seem better to account for both bar and 
arms in every disk, by measuring the mean torque over 
the whole disk: 



Qb = 



^(r)rdr 



(2) 



where rmax is the optical radius and rmin is 7% of Tmax (to 
avoid central singularities). One could also try to eliminate 
spiral structures and to preferentially focus on bars. To 
accomplish this, instead of measuring the maximal torque, 
we measure the torque in the m—2 Fourier component: if 
the potential is decomposed as 



$(r,^) = $o(r-) 



' ^mir) cos{m9 



(3) 



we define Q2{r) = 2<i>2/rli^o('')| and Q2 = maxrQ2(r). 
This eliminates spirals with more than 2 arms, thus this 
parameter may usually be regarded as a measure of bar 
strength (Combes & Sanders 1981). Finally, one may in- 
corporate the spatial extent of bars by measuring: 



Q2{r)rdr 



rdr 



(4) 



In contrast, Whyte et al. (2002), have determined the 
distribution of bar strengths from the OSU sample using 
bar shapes. They define their bar strength parameter by: 

2 r 1 1 1 

-Fbar = — arctan(&/a)^ 2 _ arctan(&/a)2 (5) 

TT L J 

where b/a is the axis ratio of the bar, i.e. the minimum 
value of the axis ratio of isophotes in an ellipse-fitting 
model (Abraham et al. 1999). Whyte et al. (2002) also 
find that nearly axisymmetric disks are rare. 

It is difficult to justify one indicator instead of another 
one; indeed it is important to test our accretion hypothesis 
with several indicators. The results are given in Table |[ 
We also show histograms of Fbar in Fig. For other in- 
dicators (Q2, Q™, Q™), histograms still show a drop at 
Q ~ 0, which is only reproduced by models with accre- 
tion. 
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Fig. 2. Histogram of bar strength Fbar- Solid: observa- 
tions, data from Whyte et al. (2002) - Dashed: simula- 
tions with gas accretion that doubles the disk mass in 10 
Gyrs - Dotted: simulations without gas accretion. 



All these additional indicators confirm the previous re- 
sults: the properties of spiral galaxies cannot be repro- 
duced without gas accretion, but require galaxies to dou- 
ble their mass in less than a Hubble time. The best accre- 
tion rate makes the mass of spiral disks double over epochs 
spanning about 10 Gyrs. The difference between the obser- 
vations and models of closed systems - or galaxies which 
accrete gas only at a marginal rate - is so flagrant that 
such models are clearly ruled out. Only large accretion 
rates reproduce the qualitative shape of the distribution 
of gravitational torques. Observations of bar axis ratios 
give the same result. Thus, it seems reasonable to con- 
clude that galaxies accrete large amounts of gas, doubling 
their mass in less than one Hubble time, even if the ex- 
act value we have derived for the accretion rate may need 
fine-tuning in future studies. 

4.3. Other processes for bar evolution 

We have identified gas accretion as being the process re- 
sponsible for bar reformation in spiral galaxies. Galaxy 
interactions may also explain the formation of bars or 
strengthen existing bars (Gerin et al. 1990). The self- 
destruction of the first bar heats the disk up and makes it 
stable. Galaxy interactions would not cool the disk down, 
a prerequisite for bar reformation and spiral arm mainte- 
nance. In the absence of accreted gas, the disk would then 
become stable because of a large velocity dispersion, thus 
it would remain unbarred while the spiral structure would 
become weaker and weaker. Gas accretion, however, keeps 
the disk cold, and it is then much more susceptible to pe- 
riodic bar reformation and spiral structure maintenance. 

Accretion is thus much more efficient compared to 
galaxy interactions in maintaining gravitational torques 
(spirals and bars) in disk galaxies. However, including the 
influence of galaxy interactions is in preparation. This will 
allow the distribution of bar strengths from our models (as 
well as the accretion rates we derive) to be refined, but in- 
teractions do not invalidate our present conclusions. 



The case of lenticulars will be treated in a separate 
paper. SOs are stellar systems without much gas. In the 
absence of gas, the dynamics of disks is different: pure stel- 
lar bars are very robust, and can endure for one Hubble 
time (Combes & Sanders 1981), contrary to bars in spiral 
galaxies. In SOs, bars are not destroyed, and no mecha- 
nism is needed to explain bar reformation. 

5. Conclusion 

The issue of whether observations support a galaxy dis- 
solving and reforming its bar has awaited the completion 
of near-infrared surveys. We are now able to show how 
robust this scenario is. Isolated galaxies (non-accreting 
systems) cannot reproduce the observed properties at all: 
they would become unbarred and spiral arms would dis- 
appear; many disks would then be nearly axisymmetric 
after a few Gyrs. On the contrary, spiral galaxies appear 
to be open systems that are still forming and continu- 
ously accreting mass today. We expect a doubling in disk 
mass every 10 billion years. The origin of the accreted gas 
has not been considered, but the most likely source could 
be the reservoirs of gas observed outside nearly all spiral 
disks (Sancisi 1983). Pfenniger et al. (1994) have even pos- 
tulated that the dark matter around spiral galaxies might 
be in the form of cold gas. Accretion rates from infalling 
dwarf satellites only represent a few percent of the accre- 
tion rate that we derive (Toth & Ostriker 1992), so that 
other sources of accretion must be invoked. 
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